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SUMMARY 

A  procedure  "based  on  the  Weiss inger  method  has  been  devised  so 
that  the  "basic  span  loading  and  associated  aerodynamic  characteristics 
can  he  rapidly  predicted  for  wings  having  arbitrary  values  of  sweep, 
aspect  ratio,  taper  ratio,  and  twist.  A  method  for  correcting  for 
the  effects  of  compressibility  is  given.  A  comparison  of  the  results 
of  this  method  with  that  of  lifting— line  and  liftin^naurface  methods 
indicates  that  the  accuracy  is  much  better  than  that  obtained  with 
lifting— line  methods  and  is  comparable  to  that  obtained  with  lifting- 
surface  methods. 

This  report,  together  with  NACA  TN  No.  1491,  allows  a  simple  and 
rapid  prediction  of  both  the  basic  and  additional  loading  character¬ 
istics  for  wings  of  arbitrary  plan  form.  The  characteristics  which 
can  be  found  for  a  given  wing  are  as  follows: 

1.  Span  load  distribution  due  to  twist  (this  report) 

2.  Span  load  distribution  due  to  angle  of  attack  (TN  No.  1491 ) 

3.  Induced  drag  (this  report  and  TN  No.  1491 ) 

4.  Angle  of  zero  lift  (this  report) 

5.  Lift-curve  slope  (TN  No.  1491) 

6.  Pitching  moment  at  zero  lift  (this  report) 

7.  Location  of  aerodynamic  center  (TN  No.  1491 ) 

It  Is  believed  these  predicted  values  are  valid  at  all  subcritical 
Mach  numbers  and  for  all  lift  coefficients  where  viscous  and  stall 
effects  are  negligible. 
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To  establish  the  effects  of  sweep,  aspect  ratio,  and  taper 
ratio  on  the  basic  loading  characteristics  produced  by  uniform 
twist,  the  method  presented  in  this  report  was  applied  to  a  few 
representative  wing  configurations  and  the  results  discussed. 


INTRODUCTION 

The  need  for  information  on  the  subsonic  characteristics  of 
swept  wings  to  supplement  the  limited  amount  of  existing  experimental 
data  has  lead  to  a  theoretical  study  of  their  characteristics.  As 
in  most  theoretical  studies,  the  wing  characteristics  have  been 
determined  using  the  span  loading  which  in  this  case  has  been  pre¬ 
dicted  by  the  method  of  Weissinger.  Since  these  characteristics  are 
dependent  upon  the  span  loading,  it  has  been  found  convenient  to 
study  the  characteristics  as  those  associated  with  additional— type 
loading  (i.e,,  the  loading  due  to  wing  angle  of  attack)  and  those 
associated  with  basic— type  loading  (i.e.,  the  loading  due  to  wing  twist 
or  effective  twist). 

In  reference  1,  the  Weissinger  method  was  applied  to  a  series 
of  wings  encompassing  the  probable  ranges  of  sweep,  aspect  ratio, 
and  taper  ratio  to  determine  the  wing  characteristics  associated 
with  additional-type  loading.  The  results  (including  span  load 
distribution,  spanwise  center  of  pressure,  lift— curve  slope,  and 
aerodynamic  center)  are  presented  in  graphical  form  as  a  function 
of  wing  plan  form. 

The  present  report  Is  an  extension  of  reference  1  to  facilitate 
determination  of  the  wing  characteristics  associated  with  basic— type 
loading  (span  load  distribution,  angle  of  zero  lift,  and  pitching 
moment)  for  a  wide  range  of  plan  forms.  Since  the  basic  loading  is  a 
function  of  twist  as  well  as  the  plan— form  variables  (sweep,  aspect 
ratio,  and  taper  ratio),  it  seemed  impractical  to  present  loading 
characteristics  for  all  possible  wing  configurations.  Therefore,  it 
was  intended  that  this  report  should  present  a  simple  procedure 
which  would  allow  prediction  of  the  basic  loading  for  the  wide  range 
of  plan  forms  investigated  in  reference  1,  and  should  present  the 
actual  basic  loading  for  a  few  representative  configurations  to 
establish  the  effects  of  the  various  geometric  parameters.  The 
results  of  this  work  then,  together  with  reference  1,  should  enable 
a  rapid  evaluation  of  the  wing  characteristics  associated  with  both 
the  basic  and  additional  types  of  loading  for  wings  having  sweep 
angles  ranging  from  -45°  to  75°,  aspect  ratios  of  1.5  to  10,  and 
taper  ratios  of  0  to  1.5. 
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SYMBOLS 

Aerodynamic  Parameters 

spanvise  loading  coefficient  for  unit  wing  lift  coefficient 
(additional— type  loading) 

spanvise  loading  coefficient  for  unit  twist  (basic-type 
loading) ,  per  degree 

spanvise  loading  coefficient  for  additional-type  loading 


spanvise  loading  coefficient  for  basic-type  loading 

/c7  c  c7  c  \ 

gross  spanvise  loading  coefficient  (  — J 

Ycav  cav  ' 


section  lift  coefficient  for  additional-type  loading 
(local  lift)a  1 

qS  J 

section  lift  coefficient  for  basic— type  loading 
f( local  lift)h  "I 

L  qs  J 


gross  section  lift  coefficient  (c7  +  c7  ) 

‘■a  LV 


wing  lift  coefficient  (  J 

V  qS  / 

rate  of  change  of  lift  coefficient  with  angle  of  attack 
measured  at  zero  lift 


dimensionless  circulation 


circulation,  feet  squared  per  second 

induced  drag  coefficient  [  induced  drag 

\  qS 
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pitching—moment  coefficient  due  to  basic  loading 
/ pitching  moment  due  to  basic  loading j 
qSc  J 

free_stream  dynamic  pressure,  pounds  per  square  foot 


Mach  number 

free— stream  velocity,  feet  per  second 


Geometric  Parameters 

angle  of  sweep  of  the  quarter-chord  line,  positive  for 
sweepback,  degrees 


aspect  ratio 


(?) 

taper  ratio 


wing  span  measured  perpendicular  to  the  plane  of  symmetry, 
feet 

wing  chord  measured  parallel  to  the  plane  of  symmetry, 
feet 


average  wing  chord 


(l)  «  feet 


mean  aerodynamic  chord 
wing  area,  square  feet 


b/2 

Jo 

W2  O  ay 

Jo 


,  feet 


(ar)  1 angle  of  attack  for  zero  lift  of  the  three-quarter-chord 

0  point  of  the  root  section  mean  line,  radians 

~  1 angle  of  attack  of  the  three-quarter-chord  point  of  the 

spanvise  station  V  section  mean  line,  radians 


(°V>0 


1the  values  of  ofy  for  zero  net  lift  on  the  wing,  radians 


3-All  angles  are  measured  in  a  plane  parallel  to  the  plane  of  symmetry. 
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Hwist  of  section  mean  line  relative  to  the  ving  root 
measured  at  the  three-quarter— chord  point  for  the 
spanwise  station  v  f ( Ov ) 0  —  (ar)0]  ?  radians  unless 
noted  otherwise 

1 twist  of  the  tip  section  mean  line  relative  to  the  wing 
root  measured  at  the  three-quarter-chord  point,  degrees 

lateral  coordinate  measured  from  wing  root  perpendicular 
to  the  plane  of  symmetry 

dimensionless  lateral  coordinate  (  —  ^ 

\b/2  / 

trigonometric  spanwise  coordinate  (cos-1r)) 

coefficient  indicating  the  influence  of  circulation  G  at 
station  n  on  the  downwash  angle  at  control  point  V 
where  the  location  of  n  is  defined  by  tj  =  cos  £2.  and 
the  location  of  V  is  defined  by  r)  =  cos 

8 


Subscripts 

parameters  associated  with  additional  loading 
parameters  associated  with  basic  loading 
equivalent  geometric  parameters 
value  of  parameter  at  zero  lift 
integers  defining  specific  span  locations 


PROCEDURE 


Development  of  Method 

In  the  Weissinger  lifting-line  method,  which  has  been  previously 
discussed  in  references  1  and  2,  and  is  used  herein,  the  gross 

1See  footnote  1,  p.  4. 
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circulation  (representing  additional  and/or  basic  loading)  is 
considered  concentrated  into  a  lifting  line  lying  along  the  quarter- 
chord  line.  The  boundary  condition,  fixing  the  spanwise  strength 
distribution  of  the  circulation,  requires  that  the  downwash  of  this 
lifting  line  and  its  system  of  trailing  vortices  produce  at  points 
along  the  three— quarter— chord  line  a  downwash  angle  equal  to  the 
slope  of  the  wing  section  mean  line  at  these  points.  For  the 
general  case,  the  boundary  condition  is  usually  applied  at  seven 
spanwise  stations  distributed  across  the  total  span  so  that  seven 
simultaneous  equations  can  be  formed  involving  unknown  gross 
loadings  at  the  seven  spanwise  stations.  However,  for  the  symmetric 
loading  case,  which  is  considered  herein  (fig.  l)  and  in  reference  1, 
the  loading  and  control  points  on  only  half  of  the  wing  need  be  con¬ 
sidered  so  that  determination  of  the  span  loading  requires  the 
solution  of  only  four  simultaneous  equations  of  the  form 


av 


4 

a^n  °n 

n=l 


(1) 


where 

av 


the  section  mean  line  angle  of  attack  at  the  three  quarter 
chord  point  for  the  spanwise  stations  V ,  where  the 
locations  of  V  are  defined  by  q  =  =  cos  -g- 


av,n 


n 


influence  coefficients  involving  the  purely  geometric  wing 
parameters  A,  A,  and  A 


the  unkown  dimensionless  circulation  ^ 
stations  n  where  the  locations  of  n 
7  nit 

=  V2=  cos  T 


at  the  spanwise 
are  defined  by 


In  this  form  the  set  of  equations  may  be  used  to  obtain  the  gross 
loading  on  any  wing  for  which  the  aspect  ratio,  sweep,  taper  ratio, 
twist,  and  angle  of  attack  are  specified. 


Past  experience  has  indicated  that  the  gross  loading  can  be 
better  studied  if  broken  down  into  the  basic  and  the  additional  type 
of  loading.  The  basic  loading  is  that  existing  with  zero  net  lift 
on  the  wing  and  is  due  to  twist  or  effective  twist  (e.g. ,  partial- 
span  flap  deflection  or  spanwise  change  in  camber)  of  the  wing  chord 
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plan®.  In  contrast,  the  additional  loading  is  that  producing  net 

®  winS  and  is,  in  effect,  the  loading  existing  on  an 
untwisted  and  uncambered  wing.  Prom  this,  it  follows  that  the 
basic  loading  is  a  function  of  the  variation  in  section  angle  of 
attack  across  the  wing  and  is  independent  of  the  wing  angle  of 
attack;  whereas  the  additional  loading  is  a  function  of  wing  angle 
of  attack  and  is  independent  of  the  variation  in  section  angle  of 
attack  across  the  wing.  ^ 

in  r-,,SndeteiTi+ing  additional  loading  characteristics  presented 
n  reference  1  then,  it  was  not  necessary  to  consider  variation  of 
the  angle  of  attack  across  the  span  (i.e.,  it  was  assumed  0^  = 

=_av=3  =  av=4)  and  H  was  possible  to  use  equation  (l)  in* the 
following  form: 


n=l 


However,  in  determining  the  basic  loading,  the  angle  of  attack 
varies  across  the  span  and  equation  (l)  must  be  used  as  given. 
Although  the  variation  of  a  across  the  span  is  known  from  wing 
geometry  for  any  wing  angle  of  attack,  its  value  for  the  unique 
condition  of  zero  net  lift  on  the  wing  is  required  and  this,  of 
course,  is  unknown.  Thus,  there  results  a  set  of  four  simultaneous 
equations  which  involve  eight  unknowns  —  the  values  of  (av)0  and 
Gn*  Since  the  twist  distribution  is  known,  three  of  these  unknowns 
can  be  eliminated  by  the  following  expression: 

(Vo  *  (ar'o  +  ev 

where  (ar)o  the  angle  of  attack  of  the  root  mean  line  at  the 
three-quarter— chord  point  for  CL  =  0,  and  ev  is  the  wing  twist 
relative  to  the  root.  This,  then,  reduces  the  number  of  unknowns 
in  the  set  of  equations  to  five  and,  hence,  a  fifth  equation  is 
required  to  permit  a  solution. 


The  fifth  equation  is  obtained  from  the  expression  for  the 
total  lift  on  the  wing.  This  expression  is  given  as  equation  (C5l) 
of  reference  2  and  where  the  seven-point  (symmetric  loading) 
solution  is  used  reduces  to 
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_  rtA 

L"T 


Gn  sin  <Pn 


(2) 


Since  cpn  =  and  since  for  basic  loading  CL  =  0,  the  only 

uhkowns  in  equation  (2)  St'ionf^at 

appear  in  equation  U)«  ’  and  G4  .  In  the  expanded 

contain  the  unknowns  (ar)0#  ^2  >  «3 '  4 

form  the  equations  are: 

s€l  =  -  (ar)0  +  a  &  +  a1>2  G2  +  a1>3G3  +  a1)4G4 

«8  -  -  M0  *  %1Q1  +  a2,2  G2  +  a2,3G3  +  H,*** 

e3  =  -  (or)o  +  a3  Gx  +  ag^2  G2  +  &3^3G3  +  %4  4 

0  -  -  (ar)0  +  a4<1Gx  +  a4>2  G2  +  a4)3G3  +  a4j4  4 

Cl=0  =  O.765G1  +1.414G2+  1.848g3  +  G4 

Solution  of  this  sot  of  equations  vill  gi-  * ho  angle  of  attach 

of  tha  mean  line  ^^"^afions  ,  -  0.92U, 

and  the  loads  Gi,  G2,  ^4 

0.707,  0.383,  and  0. 

Application  of  Method 

naaic  loading  and  angl^ggU^- 

h“ S'deXmina  is  loading  on  a  ving  having  a 


equations  should  he  written  tan[e1+(or)0]  -  ai,iGx  x’f. 

Early  in  the  d«riTrti«m 1  Us  ^66^1^0  substitute 

slopes  were  small  ^  £  avoid  serious  error  in 

the  angle  iT1  **a  ^  twi6t  is  large,  the  equations  should  he 

those  cases  where  the  twist  is  larg  ,  *  radians.  If  it 

HBI 

should  he  used  for  angles  over  20  . _ _ _ _ _ 


MCA  TN  No.  1772 


9 


given  twist  distribution  or,  conversely,  to  determine  the  twist 
distribution  to  provide  a  given  basic  loading.  In  either  case  the 
solution  of  the  five  simultaneous  equations  is  a  relatively  simple 
matter.  The  most  time-consuming  and  laborious  portion  of  the  process 
is  that  of  computing  the  values  of  the  av  n  coefficients  from  the 
geometry  of  the  wing.  This  process  is  fully  outlined  in  reference  2. 
Since  these  coefficients  are  a  function  of  wing  geometry  alone, 
however,  it  is  only  necessary  to  compute  them  once  for  each  plan 
form  to  study  any  variation  of  twist  or  camber  on  the  plan  form. 

These  computations  have  been  made  for  the  range  of  plan  forms  shown 
in  figure  2  and  are  presented  in  table  I.  (it  should  be  noted  that 
figure  2  shows  only  the  range  and  not  the  total  number  of  plan  forms; 
approximately  200  wings  were  considered  altogether.) 

Thus,  with  the  aid  of  table  I,  the  problem  of  determining  the 
basic  loading  and  the  angle  of  zero  lift  for  any  plan  form  given  in 
table  I  is  reduced  to  the  following  simple  steps: 

1.  Insertion  (in  equation  (3))  of  the  given  values  of  twist 

for  the  four  spanwise  stations 

2.  Insertion  of  the  values  of  av^n  obtained  from  table  I 

3.  Simultaneous  solution  of  the  five  equations 

The  resulting  loading  coefficients  may  be  put  in  the  more  convenient 

c7  c 
t-b 

coefficient  form  — -  by  the  following  conversion: 

cav 


c  _  g 

To  aid  in  fairing  the  loading  curve,  values  of  -=&—  at  intermediate 

cav 

span  stations  may  be  obtained  through  use  of  the  interpolation 
function  (equation  (A 6))  given  in  reference  1.  Determination  of  the 
twist  for  a  given  loading  is,  of  course,  a  simple  inversion  of  this 
process. 

As  will  be  evident  later,  it  is  not  necessary  to  use  the  exact 
plan  form  to  obtain  a  good  approximation  of  the  basic  loading  on  a 
given  wing.  Consequently,  if  the  basic  loading  is  desired  on  a  plan 
form  between  those  given  in  table  I,  it  is  generally  acceptable  to 
use  the  coefficients  av^n  given  in  table  I  for  a  plan  form  which 
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most  closely  approximates  the  desired  plan  form.  If  a  more  refined 
estimate  of  loading  is  required,  experience  has  indicated  that  it  is 
better  to  determine  the  loading  for  the  bracketing  plan  forms  and 
interpolate  the  loading  rather  than  to  determine  the  loading  for  the 
interpolated  values  of  the  coefficients  Qy.n 

Consideration  must  be  given  to  the  number  of  significant  figures 
retained  in  the  solution  of  the  five  simultaneous  equations. 

Actually  the  number  of  significant  figures  required  depends  to  a 
large  extent  on  whether  or  not  the  equations  are  ill— conditioned} 
however,  it  is  noteworthy  that  in  solving  for  the  basic  loadings 
presented  in  this  report,  none  of  the  sets  of  equations  appeared  to 
be  ill-conditioned.  Beyond  possible  effects  of  ill-conditioning, 
however,  to  be  strictly  correct,  the  number  of  significant  figures 
to  be  retained  should  be  examined  at  each  step  of  the  computations 
to  maintain  a  given  accuracy.  When  computations  are  made  either 
longhand  or  with  a  slide  rule,  this  procedure  can  be  followed,  but 
when  computing  machines  are  used  it  becomes  more  practical,  even 
though  not  rigorously  correct,  to  maintain  a  given  number  of  decimal 
places  throughout  the  computations .  In  an  effort  to  establish  the 
number  of  decimal  places  required,  a  number  of  computations  were 
made  using  five  places  and  then  four,  three,  etc.,  for  each  of 
several  plan  forms.  The  results  obtained  were  compared  and  it  was 
concluded  that  satisfactory  accuracy  could  be  had  if  the  value  of  ev 
in  radians  were  given  to  four  decimal  places,  if  the  coefficients 
ay  n  were  tabulated  and  used  to  two  decimal  places,  and,  if  in  the 
solution  of  simultaneous  equations,  five  decimal  places  were  retained. 

Local  lift.  Induced  drag,  and  pitching  moments.-  With  the  basic 

Ci  c 

loading  coefficient  — —  as  determined  above  and  the  additional 

cav 

loading  coefficient  tl&L  as  determined  by  method  given  in  reference  1, 

cav 

other  wing  characteristics,  such  as  section  gross  local  lift  coeffi¬ 
cient  cz,  induced  drag  coefficient  Cj)±,  and  wing  pitching-moment 
coefficients  Cm  are  easily  obtained.  The  gross  value  of  c j  at 
any  angl e  of  attack  is  determined  as  follows: 

cl  =  cZa  +  clb  '  (5) 

where  c7  is  the  lift  coefficient  due  to  additional-type  loading 
from  reference  1  and  c7  is  the  lift  coefficient  due  to  the  basic 
type  loading  as  determined  by  the  relation  # 
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(i+a) 

2  [1hi(1-*)] 


In  contrast  to  lift  coefficient  and  loading,  the  induced  drag 
cannot  be  determined  for  additional  and  basic  loading  separately  and 
then  be  summed  to  get  the  total  induced  drag,  rather  the  induced 
drag  must  be  determined  from  the  total  loading  distribution. 

Equation  (3)  of  reference  1  has  been  modified  to  give  the  induced 
drag  for  the  general  case.  The  induced  drag  coefficient  is  then 
given  by 


°Di  =  [  *i2  +  ^  +  -f 


k4  (0.0561  kx 


+  0.7887  k  )  -  k2  (0.7352  kx  +  0.8445  k3)] 


where 


kn  — 


The  ability  of  the  Weissinger  method  to  enable  good  predictions 
of  the  effect  of  load  on  wing  pitching  moment  has  been  shown 
(reference  l)  to  be  the  result  of: 

1.  The  accurate  prediction  of  spanwise  distribution  of  load 

2.  The  predominant  effect  that  spanwise  distribution  of  load 

has  on  pitching  moments  of  swept  wings  as  compared  to  the 
effect  of  chordwise  distribution  of  load 

Thus,  even  though  in  the  Weissinger  method  the  basic  load  distribu¬ 
tion  is  concentrated  along  the  quarter-chord  line,  the  method  should 
allow  good  predictions  of  the  effect  of  basic  loading  on  the  pitching 
moment  of  swept  wings. 


The  expression  for  due  to  basic  loading  has  been  derived 

from  equation  (A*0  of  reference  1.  Thus, 

Cntb  =  tan  A  (0.1384  klfe  +  0.1975  k2b  +  0.1351  k3b  +  0.0159  k4b)  (8) 


c 
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It  should  be  noted  that,  since  the  pitching  moment  due  to  basic 
loading  is  the  result  of  a  loading  couple,  the  value  of  the  pitching 
moment  due  to  basic  loading  is  independent  of  the  location  of  the 
moment  reference  center.  To  obtain  the  gross  pitching-moment  coeffi¬ 
cient  for  a  wing,  the  pitching  moment  due  to  additional— type  loading 
(reference  l)  may  be  added  directly  to  the  pitching  moment  due  to 
twist  given  by  equation  (8). 


Effects  of  compressibility.—  A  means  of  correcting  the  wing 
characteristics  associated  with  additional— type  loading  for  compress¬ 
ibility  was  given  in  reference  1.  This  essentially  consists  of 
translating  the  effect  of  compressibility  into  an  effective  change 
in  plan  form  in  addition  to  the  well-known  increase  in  section 
pressures.  These  principles  should  apply  equally  well  to  basic 
loading  characteristics.  However,  it  should  be  noted  that  in  the 

czc 

case  of  additional  loading  the  loading  coefficient  -  was  not 

CLcav 

a  function  of  angle  of  attack  and  consequently  the  effects  of  Mach 
number  on  loading  were  shown  only  as  changes  in  load  distribution 


and  not  as  changes  in  the  average  value  of 


C-j  c 


The  change  in 


^Lcav 

average  value  of  loading  was  in  effect  absorbed  in  changes  in  Ct 

c2v>c 

In  contrast,  the  basic  loading  coefficient  — —  is  definitely  a 

cav 


function  of  the  local  angle  of  attack  (twist)  and  is 
function  of  Mach  number  just  as  is  lift-curve  slope. 

°l  0 

obtain  the  value  of  — in  compressible  flow,  it 

cav 


,  therefore,  a 
Therefore,  to 

is  necessary  to: 


1.  Determine  the  value  of 


c-,  c 


av 


for  the  given  twist  and  the 

equivalent  plan  form  given  by  \  Ae  =  A  ,/l-M2 

tan  A 


and  tan  Aq  = 


2.  Multiply  value  of 


ci-bc 


cav 


obtained  by 


Vi"**2 
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RESULTS  AND  DISCUSSION 
Evaluation  of  Method 

There  is  a  scarcity  of  experimental  basic  loading  data  and 

^7+?Valrti0n  °f  the  accurac7  of  the  Weissinger 
method  in  predicting  basic  loading  must  be  indirect.  Bothrefer- 

ences  1  and  2  proved  the  Weissinger  method  to  be  very  accurate  L 

ther^fn1118  additi0ml~t7pe  loadlng,  and  similar  accuracy  should 
therefore ,  be  expected  with  regard  to  the  basic  loading. 

givenTinafi^rIT?f  eval^tion  of  the  ^thod,  a  comparison  is 
L®. in  figure  3  between  the  basic  loadings  obtained  by  the 

inger  method,  the  method  of  reference  3  and  the  Falkner  method 

0^0??.^ LS*?  hSVin?  an  SSpeCt  rati°  °f  6-°  a  *£r^S£ 

. ,  rom  reference  3  were  used  for  comparison  since 

hey  are  well  known  and  have  been  widely  used.  The  Falkner  method 
h?+USed  because  it  is  a  lifting-surface  method  and  should  give 

^6r  ^Sn  6lther  the  miM  ot  reference  3  (a  liftiL- 

line  method)  or  the  Weissinger  method  (a  modified  lifting-line 

”®^°d)*  pr®sented  111  figure  3,  the  loading  obtained  by  refer- 
nce  3  is  in  serious  disagreement;  whereas  the  Weissinger  loading 
shows  reiativeiy  good  agreement  with  that  obtained  by  the  Falkne? 

factsf  *  rS8Ulta  "&***>»*  «  the  basis  It  the  filing 


.  C?'n  be  readily  shown  that,  even  on  high— aspect— ratio 

wings,  the  introduction  of  twist  results  in  large  induction  effects, 

!her?  induction  effects  are  large,  as  for  example  the 
3  f  ^dacti°n  °n  the  lift-curve  slope  of  low-aspect-ratio 
wings,  it  has  been  often  demonstrated  that  unmodified  lifting-line 
theory  will  not  yield  accurate  results.  lifting-line 

,  ,  3*  In  reference  1  it  was  shown  that  the  Weissinger  method 

theCun^odifT°di+'he<a  lifting“line  method^  overcomes  the  weakness’ of 
the°ry  and  7lelds  results  on  low-aspect— ratio  wings 
comparable  m  accuracy  to  that  obtained  with  lifting-surface  theory. 

th*  the^ foregoing  comparisons,  it  is  believed  that  (l) 

the  basic  loading  characteristics  of  unswept  wings  can  be  predicted 

us  w^e  ?+r  aCra7  by  USlnS  the  method  than  by 

using  the  results  of  reference  3,  and  (2)  that  the  Weissinger  method 

wingawithes0-Ff*?r?di+tin8  thS  basic  loading  characteristics  on  any 
wing  with  sufficient  accuracy  for  preliminary  design  analysis. 
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Effect  of  Plan-Form  Variation  on  the 
Basic  Loading  Characteristics 

To  study  the  effects  of  plan-form  variation  on  the  basic 
loading  characteristics  of  uniformly  twisted  wings,  the  character¬ 
istics  of  a  representative  group  of  wings  (see  shaded  wings,  fig.  2) 
having  unit3  washout  have  been  computed  and  are  presented  in  figures 
4  to  9.  The  basic  loading  characteristics  considered  are  the 
c  c 

loading  —  (figs.  4  to  7) ,  the  pitching  moment  due  to  twist  Cm 
cav  .  . 

(fig.  8),  and  the  angle  of  attack  at  the  root  for  zero  lift  (or)Q 

(fig.  9). 


Magnitude  and  spanwise  distribution  of  load.-  Examination  of 
figures  4  to  7  reveals  that  the  aspect  ratio  influences  only  the 
magnitude  and  is  in  fact  the  predominate  influence  on  the  magnitude. 
Reductions  in  aspect  ratio  from  6.0  to  3«5  and  1.5  result  in  approxi¬ 
mately  35-percent  and  70-percent  reductions,  respectively,  in  load 
due  to  twist  for  either  the  unswept  or  45°  swept-back  wings  (fig.  5). 

Sweep,  either  forward  or  back,  tends  to  reduce  the  magnitude  of 
loading,  although  appreciable  reductions  are  produced  only  by  sweep 
angles  greater  than  45°  (fig.  4) .  Sweep  also  affects  the  load 
distribution  such  that  the  load  on  the  outer  section  of  the  wing  is 
shifted  inboard’ by  sweepforward  and  toward  the  tip  by  sweepback; 
this  is  similar  to  the  effect  of  sweep  on  the  additional-type 
loading.  Since  increase  in  aspect  ratio  magnifies  the  loading,  it 
also  magnifies  the  effects  of  sweep  on  the  loading  as  is  shown  m 

figure  4. 

As  shown  in  figures  6  and  7,  taper  ratio  has  little  effect  on 
the  magnitude  of  basic  loading;  and  variations  in  taper  ratio,  for 
taper  ratios  larger  than  0.5,  have  little  effect  on  the  load 
distribution.  However,  for  taper  ratios  less  than  0.5>  tbe  loading 
on  the  outer  section  of  the  wing  shifts  inboard,  These  effects  of 
taper  ratio  on  loading  are  magnified  by  increases  in  aspect  ratio. 


3 In  this  case,  1°  was  chosen,^ and  for  any  larger  amount  of  twist 
the  effects  are  proportional  within  the  limits  of  footnote  2, 
page  8.  _ __ _ _ _ 
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Pitching  moment; .  —  That  the  pitching  moment  due  to  twist  is 
primarily  a  function  of  sweep  and  aspect  ratio  is  shown  in  figure  8. 
The  magnitude  of  the  pitching  moment  increases  as  either  aspect 
ratio  or  sweep  is  increased  so  that  pitching— moment  coefficients  as 
large  as  0.008  for  1®  of  twist  exist  on  wings  having  large  aspect 
ratios  and  sweep  angles.  The  effect  of  taper  ratio  is  relatively 
small,  the  greatest  being  evidenced  at  the  small  values  of  taper 
ratio.  For  example,  reducing  the  taper  ratio  from  0.5  to  0  reduces 
the  pitching  moment  due  to  twist  about  30  percent. 


Angle  of  zero  lift.-  Although  the  effects  of  plan  form  on  the 
angle  of  zero  lift  (or)0  may  not  be  very  important,  some  of  the 
trends  indicated  in  figure  9  are  of  interest.  For  the  range  of  plan 
forms  represented  in  figure  9,  the  angle  of  zero  lift  did  not  vary 
more  than  27  percent.  This  is  small  compared  to  the  effects  of  plan 
form  on  the  magnitude  and  distribution  of  loading  and  on  the  pitch¬ 
ing  moment.  In  contrast  to  the  small  effect  of  taper  noted  previ¬ 
ously,  taper  ratio  appears  to  be  the  predominant  influence  on  (c^)  , 

particularly  at  large  aspect  ratios  and  large  sweepback.  The  effect 
of  aspect  ratio  and  sweep  are  secondary  but  not  negligible. 


Consideration  of  Twist  in  Swept-Wing  Design 

The  methods  of  this  report  enable  a  detailed  study  of  two  seri¬ 
ous  problems  associated  with  the  use  of  swept  wings:  First,  the  use 
of  twist  to  control  section  lift  distribution  and  hence  stall;  and, 
second,  the  degree  to  which  bending,  since  it  introduces  twist, 
affects  the  aerodynamic  characteristics.  These  two  problems  and  the 
application  of  the  methods  of  this  report  to  this  study  are  discussed 
in  the  sections  immediately  following. 

Twist  for  separation  control.-  The  induction  effects  on  swept 
wings  are  such  that  large  angles  of  attack  are  induced  near  the  root 
of  swept-forward  wings  and  near  the  tip  of  swept-back  wings.  Conse¬ 
quently,  for  untwisted  wings  the  sections  at  the  root  and  tip, 
respectively,  reach  their  cimny  before  the  rest  of  the  wing  and  at 
a  relatively  low  angle  of  attack  of  the  wing.  The  resulting  local 
separation  produces  the  poor  characteristics  (large  drag  rise  and 
large  fore-and-aft  movement  of  the  aerodynamic  center  which  occur 
at  relatively  low  lift  coefficients)  which  are  typical  of  highly 
swepo  wings.  To  remedy  these  poor  characteristics,  some  means  must 
be  provided  which  will  cause  the  flow  over  the  wing  to  separate  more 
uniformly.  This  should  be  achieved  if  the  wing  is  twisted  and/or 
cambered  so  that  all  sections  reach  their  c  at  nearly  the  same 
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angle  of  attack  of  the  wing.  One  possible  approach  would  be  to  twist 
the  wing  to  provide  nearly  uniform  spanwise  distribution  of  4c j  and, 
hence,  uniform  separation,  and  to  camber  all  sections  to  increase 
cLmax- 

An  indication  of  the  amount  of  twist  required  to  produce  uniform 
distribution  of  cj  for  a  given  can  be  had  by  (l)  using  refer¬ 

ence  1,  determine  the  Cj  distribution  for  the  given  CL  (additional- 
type  loading),  (2)  using  this  Cj  distribution,  determine  the  basic 
loading  required  to  give  uniform  cj  distribution,  and  (3)  using  the 
method  of  this  report,  compute  the  twist  distribution  required  to 
produce  this  basic  loading. 

For  purposes  of  illustration,  the  foregoing  procedure  has  been 
applied  to  a  wing  having  45°  of  sweepback,  an  aspect  ratio  of  6.0, 
and  a  taper  ratio  of  0.5*  (See  fig.  10.)  It  can  be  seen  that  the 
additional-type  loading  produces  relatively  large  values  of  cj 
over  the  outboard  sections.  The  basic  loading  required  to  give  uni¬ 
form  c l  distribution  of  1.0  is  shown  shaded,  and  the  basic  loading 
ordinates  to  be  used  in  equation  (3)  in  order  to  determine  the 

required  twist  are  labeled  (cj,)  >  ^cl>)  *  e"kc*  The  twist,  deter— 

4  9 

mined  from  solution  of  equation  (3)  to  satisfy  the  condition  of  uni¬ 
form  load,  is  shown  on  the  lower  half  of  the  figure.  Such  a  vari¬ 
ation  would  be  difficult  to  build,  and  hence  in  the  practical  appli¬ 
cation  some  compromise  twist  variation  would  probably  be  chosen. 

The  effect  of  compromising  to  the  extent  of  using  a  linear  variation 
in  twist  is  also  shown  in  figure  10.  Although  a  uniform  cj  dis¬ 
tribution  is  not  provided,  the  distribution  is  such  that  the  char¬ 
acteristics  of  the  wing  at  higher  lift  coefficients  should  be  signi¬ 
ficantly  improved.  Undoubtedly  the  amount  of  twist  and  camber  would 
also  be  compromised  to  some  extent  to  provide  good  characteristics 
at  high  Mach  numbers;  however,  it  is  likely  that  relatively  large 
amounts  of  camber  and  twist  could  be  tolerated  on  highly  swept  wings 
before  the  characteristics  at  high  Mach  numbers  were  jeopardized. 


4It  is  recognized  that  the  optimum  spanwise  distribution  of  cz 
may  not  be  a  uniform  distribution  and  that  in  the  practical 
application  better  stalling  characteristics  will  probably  be 
exhibited  if  the  c,  at  the  midsemispan  is  somewhat  greater 
than  at  either  the  root  or  tip  sections.  It  is  apparent  that 
further  experimental  data  are  needed  to  establish  the  optimum 
Cj  distribution  for  swept  wings. 
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Twist  due  to  wljng  deflection.—  Deflection  of  wings,  whether 
unswept  or  swept,  can  he  considered  as  composed  of  tending  and  tor¬ 
sional  components.  Torsional  deflection  of  either  unswept  or  swept 
wings  obviously  produces  twist  and,  hence,  affects  the  basic  loading. 
On  unswept  wings,  bending  produces  only  an  Increment  In  the  dihedral 
angle;  whereas  on  swept  wings,  bending  produces  changes  In  section 
angle  of  attack  (effective  twist)  as  well  as  dihedral.  Whether  or 
not  the  twists  due  to  bending  and  torsional  deflections  are  additive 
or  canceling  depends  upon  the  direction  of  sweep  and  the  character¬ 
istics  of  the  wing  structure.  For  a  given  wing,  the  magnitude  of  the 
twist  due  to  wing  deflection  is  a  function  of  the  gross  loading  as 
well  as  the  structural  stiffness  and,  therefore,  varies  with  accel¬ 
eration  in  gusts  or  maneuvering  flight.  Thus,  the  aerodynamic  char¬ 
acteristics  of  the  wing  in  maneuvering  flight  are  likely  to  differ 
greatly  from  those  of  the  wing  in  steady  flight. 

For  swept  wings,  the  pitching  moment  resulting  from  this  twist 
directly  affects  the  trim  and/or  stability  of  the  airplane  and  should 
therefore,  be  given  serious  consideration.  A  qualitative  analysis 
has  shown  that  for  a  flexible  wing  in  steady  flight,  the  pitching- 
moment  increments  due  to  flexure  result  primarily  in  changes  in  trim 
but  may  also  cause  decrements  in  stability  which  increase  in  magni¬ 
tude  as  lift  coefficient  is  decreased.  In  accelerated  flight,  the 
flexible  wing  will  probably  experience  serious  decrements  in  stabil¬ 
ity  at  all  lift  coefficients. 

These  changes  in  trim  and  stability  can  be  evaluated  through 
use  of  equations  (3)  and  (8)  once  the  structural  stiffness  of  the 
wing  in  both  bending  and  torsion  (hence,  the  twist  components)  has 
been  determined.  As  can  be  seen  in  figure  8,  the  effects  of  twist 
on  trim  and  stability  are  very  dependent  on  the  sweep  and  aspect 
ratio  of  the  wing,  and  to  same  extent,  dependent  on  the  taper  ratio, 
Beducing  either  aspect  ratio  or  sweep  will,  therefore,  minimize  the 
influence  of  twist  and  wing  bending  on  the  trim  and  stability  of  the 
airplane. 


CONCLUDING  REMARKS 

A  method  for  predicting  effects  of  twist  on  the  span  loading 
and  associated  characteristics  for  a  wide  range  of  plan  forms  has 
been  presented.  Comparison  of  the  loadings  obtained  by  this  method 
and  those  obtained  by  lifting-surface  (Falkner)  and  lifting-line 
(NACA  TR  No.  572)  methods  indicates  that  use  of  this  method  results 
in  accuracy  much  better  than  that  obtained  with  lifting-line  methods 
and  accuracy  comparable  to  that  obtained  with  lifting-surface  methods 
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It  is  believed,  therefore,  that  (l)  this  method  is  capable  of  pre¬ 
dicting  the  basic  loading  characteristics  on  any  wing  with  sufficient 
accuracy  for  preliminary  design  analysis,  and  (2)  the  method  should  be 
particularly  useful  in  determining  the  twist  required  for  separation 
control  and  in  evaluating  effects  of  aeroelastic  deformation. 

To  establish  the  effects  of  sweep,  aspect  ratio,  and  taper  ratio 
on  the  basic  loading  characteristics  produced  by  uniform  twist,  the 
method  has  been  applied  to  a  few  representative  wing  plan  forms. 

Same  of  the  trends  noted  were: 

1.  The  magnitude  of  basic  loading  due  to  uniform,  twist  is 
primarily  a  function  of  aspect  ratio;  however,  the  influence  of  sweep 
also  becomes  important  for  sweep  angles  beyond  45°.  Taper  ratio  has 
little  effect. 

2.  The  pitching  moment  due  to  uniform  twist  is  a  function  of 
both  sweep  and  aspect  ratio  and  is  likewise  little  affected  by  taper 
ratio. 


3.  The  of  zero  lift  of  a  uniformly  twisted  wing  is  a 

function  of  taper  ratio  as  well  as  sweep  and  aspect  ratio. 


Ames  Aeronautical  Laboratory, 

national  Advisory  Committee  for  Aeronautics, 
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Note:  Basic  hading  characteristics  of  shaded  wings  are  discussed  in  text. 
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Figure  2—  Approximate  range  of  ptanforms  for  which  values 
of  coefficients  0  ore  presented  in  table  I. 
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Spanwise  loading  coefficient  for  unit  twist,  ,  per  deg 
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Figure  4-  The  effect  of  sweep  on  the  basic  loading  of 
wings  having  a  taper  ratio  of  0.5 . 
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Spanwise  station,  ij,  fraction  of  semispan  ''  .Cnaca'J^’ 


Figure  5.  -  The  effect  of  aspect  ratio  on  the  basic  load¬ 
ing  of  wings  having  a  taper  ratio  of  0.5 . 


S panwise  loading  coefficient  for  unit  twist \  —  .  per  deg 

*  Cav 


MCA  TN  Wo.  1772 


;  0 


-.004 

-.008 

-.012 


A--  .5 
\  =  1.5 
\  =  0 


A-0 


O  ./  .2  .3  .4  .5  .6  .7  .8  .9 


-5 

ru:,: 


A  =45 


012 1 _ 1 _ 1 _ 1 _ 1 _ 1 _ 1 _ 1  1  1 

O  ./  .2  .3  .4  .5  .6  .7  .8  .9 

Spanwise  station,  y,  fraction  of  semispan  «=__ 


Figure  6- The  effect  of  taper  ratio  on  the  basic  loading 
of  wings  having  an  aspect  ratio  of  3.5 . 
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Figure  7rThe  effect  of  taper  ratio  on  the  basic  loading 
of  wings  having  an  aspect  ratio  of  6.0. 


Figure  8—  The  effect  of  plan  form  on  the  pitching  moment  due 
to  twist. 


Spanwj  se  station,  -q,  fraction  of  semi  span 


Figure  10. — Example  illustrating  use  of  twist  to  control  spanwise 
distribution  of  local  lift  coefficient  for  a  wing  having  45° 
of  sweepback ,  an  aspect  ratio  of  6.0,  and  a  taper  ratio  of  0.5. 


J 


p 

..  g 

PH  P 
P  to 
«p  O  P 
O  ©  P 
Ph  CD 
(Drop 


t5 

P 

t£) 

cd  o  • 

p 

P 

P 

cd  p  d 

T* 

© 

$5 

a 

cd 

•>  P  P  © 

ol 

P 

Ph 

G 

P 

Ph  cd  cd  ra 

P 

m 

P 

©  p  p  ra 

.p 

© 

©  o  P  P 

ra 

P 

o 

PD  d 

>  _  ©  o 

_Q 

O 

p 

P 

p 

ra  PD  ra  m 

5 

Ph  ra 

P 

p  ©  p 

1  p 

cd 

P 

P  p  P  d 

cd  o 


Ph 


© 


d 

© 

o 

o 

p 

Ph 


d 
©  P 
P  P 

P  H 

© 

ra  o 
©  P 
P  © 
Pi  n 


© 
P 
Pi  O 
©  © 
P  d  P 
POP 
P  © 
d  o 
P  P  © 
cd  d  P 
©  P 
*  P 

P  Pi  P 
CD 
►>a  P 


Pi 
©  ra 
P  P 


o 


>  p 
©  ra 
P  © 
P 
cd 

© 

o  p 
p  p 


d  d 
©  P 


P  p 
d  cd 
ed  P 
©  m 
P  © 


i — l 

P< 

Ph  P 
cd  © 


©  > 
P  P 


P  p  O 
O  d  cd  p 
P  cd  P  P 
©  ©  ra  cd 
tsl  P  ©  P 


g 


P 

P 

P  © 

© 

d  P 

O  P 

O  Ph 

O  O 

ra 

^  cd 

P  P 

p 

©  p 

p 

P  P 

H  05 

©  S 

d 

PO  O 

•  d 

£  3 

O 

P 

ra  p 

P 

cd  P 

P  ad 

©  w 

P 

ra 

p  P 

© 

©  p 

O  *vp  P 

P 

P  > 

P  o 

i> 

p  P 

p 

*'  cd 

P 

p  p 

m  P 

© 

^  o 

cd  cd 

a 

GO 

tt)  p 

H  P 

p  m 

p 

P 

Ph 

P  t}D 

p 

p  © 

P 

m  p 

m 

d  P 

P  O 

p  p 

ra 

cd  d 

o  © 

P  > 

p 

o 

Ph 

© 

© 

p  p 

©  ra 

p  © 

p 

w)  <d 

o  t> 

• 

p  © 

P 

cd  p 

d 

© 

cd  P 

cd 

P  P 

© 

p 

Ph  O 

P  Ph 

cd  cd 

ra 

p 

m  p 

© 

P  P 

ra 

©  © 

o  p 

3 

P 

O 

d  > 

© 

o 

o 

P  R 

p  ra 

ra 

m 

Ph 

ra  p 

> 

P  © 

p 

p 

id  P 

P 

p  P 

d 

p  o 

cd  O 

d  Ph 

d 

p 

td  © 

ra 

© 

©  p 

P  ra 

o  P 

p 

CD 

P  Ph 

O  P 

i — l 

r— i 

cd 

P 

<P  o 

2 

P 

© 

© 

o  © 

ra 

P  P 

d  P 

p  P 

© 

© 

cd  P 

©  P 

ra 

P 

£ 

d 

© 

o 

o 

p 

Ph 


© 


P 
d  o 
p  p  © 
cd  d  p 
©  P 
P 

P  PH  P 
P  © 

p  p 

rl  ri  rl 

P  p 

O  d  cd 
P  cd  P 
©  ©  ra 


td  cd 
P 


©  P 
p 


© 
o  p 


P  d 
© 


Ph  N  p  ©  p 


P 
cd  ra 


Ph  ~ 
Ph  P 
td  ra 


•H  3 


N 


